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Abstract 

In this paper we analytically estimate the magnetic field scale of 
planets with physical core conditions similar to that of Earth from a 
statistical point of view. We evaluate the magnetic field on the basis 
of the physical parameters of the center of the planet, such as den- 
sity, temperature, and core size. We look at the contribution of the 
Peltier-Seebeck effect on the magnetic field, showing that an electri- 
cal thermal current can exist in a rotating fluid sphere. Finally, we 
apply our calculations to Earth and Jupiter. In each case we show 
that the thermal generation of currents leads to a magnetic field scale 
comparable to the observed fields of the two planets. 

1 Introduction 

The Earth's magnetic field is a fascinating problem that has been faced 
by many authors and the reader can find rich literature on the topic. 
For many years the intuitive idea that the magnetic field is generated 
by heavy fluid in the center of Earth subjected to the rotational mo- 
tion of our planet, has been conjectured. Many numerical works have 
started to shed light on the possible mechanism of the generation of 
Earth's magnetic field. The basic model for the generation of Earth's 
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magnetic field or of other planets, is based upon the dynamo effect of 
a turbulent convection in rotating fluids. This idea has received much 
attention in the past few years and many numerical studies based on 
the dynamo model attempted to reproduce some of the main proper- 
ties of the magnetism of celestial bodies [1, 2, 3, 4], among them the 
phenomenon of magnetic field reversal. Magnetic field reversal, the 
phenomenon for which the positions of magnetic north and magnetic 
south are interchanged, is another important feature of the terrestrial 
magnetic field that has been studied intensively, and recently a similar 
phenomenon has been reproduced in the laboratory [5] . 

In a previous work [6] the authors showed that a magnetic field can 
be generated in the laminar region of a fluid velocity under the condi- 
tion p = rjafi, but such a condition is far from the usual condition of 
the celestial body and, more important, its magnitude would be of the 
order of B ~ riRy/pjlo where O is the rotational velocity, R is the ra- 
dius of the outer core, and p is the density of the fluid in the outer core 
of Earth [7]. Inserting Earth's parameters would imply an intensity 
field of S ~ 10^ gauss that is very far from Earth's actual magnetic 
field value. Independent from the path that the above value has been 
obtained, it is a fact that the magnetic field scale B ~ QR^ppo can 
be deduced by magnctohydrodynamic equations. This is an indica- 
tion that the Earth's magnetic field could have a different origin that 
is not strictly dynamic, since the value appears to be too high when 
compared with the observed field. The above discussion leads us to 
the following question: Does there exist for Earth a magnetic field 
scale as a function of the physical system parameters? We shall show 
that there exists such a characteristic magnetic field with an intensity 
that is very close to the actual value of Earth's field. 

2 Magnetohydrodynamic equations 

Let us consider the set of equations for a plasma with finite conduc- 
tivity and constant density [8, 9] 



-v + (v.V)v 



= -VP + [VxH] xB-hf + 0- (1) 
V-v = (2) 
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^ = V X [v X B] + — V^B (3) 
ot fia 

V-B = 0, (4) 

where v is the flow velocity, H is the magnetic field related to the 
magnetic induction B via the relation B = /xH, P is the pressure of 

the gas, and p is the mass density. The gravity force density, f , takes 
the form f = p'V^ . The vector cr is defined through its components 
as follows [10] 

< = (5) 



dxk ' * \dxk dxi 

where is the viscous stress tensor, and rj is the coefficient of vis- 
cosity which is assumed constant. We also used the convention of 
dropping the symbol of sum for the repeated indexes. The current 
density J is given by the constitutive relation J = (7(E + v x B) 
where a is the electrical conductivity of the fiuid. 

The dynamic system has to be implemented using the equation of 
heat transfer in magnetohydrodynamics [8] 

where Cp is the specific heat at constant pressure, k is the thermal 
conductivity, T is the temperature of the fluid, and Q is the quantity of 
heat generated by external sources of heat contained in a unit volume 
of the fluid per unit time. 

Let us deflne a dimensionless velocity u = v/(nR), a dimension- 
less operator V' = i?,V, and a dimensionless time r = Vlt. Combining 
(1) and (5) we can rewrite the hydromagnetics equations as 



^u + (u.V)u = -^Vb + ^) + J-[Vxb]xb + ^V^(^) 

V-u = (8) 

* = Vx[uxbl + ^V^b (9) 

V-b = (10) 

where we deflned the dimensionless quantities p = P/Pq, iJj = pip/Po, b 
B/Bq with Pq as a characteristic pressure, and Bq as a characteristic 
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induction field. For brevity we redefined V' as V. Finally we defined 
the numbers 



Re = 



Rp 



Po 



Rb = A* 



pn^R^ 



Rm = crjJL^R 



,2 



V 



where Re is the Reynolds number, Rb the magnetic force number, 
and Rm = (t/ji^R^ is the Reynolds magnetic number. 

As many authors have pointed out the analytical solution of this 
equation is a very hard task and only a few exact or approximate 
cases are known (see for example [11] for a review). Our main aim is 
to find an estimation for Bq as a function of the physical parameters 
of the center of Earth (such as density, temperature, etc.) without 
necessarily solving the system (7)- (10). 



3 Remarks on field velocity components 
for a rotating sphere 



In this section we will show that for a rotating sphere, in general the 
component of the velocity is not enough to describe a time dependent 
motion, with it being understood that = QrsinO, where Q is the 
angular velocity, is an exact solution. Note that in Ref. [12] is shown 
evidence that the inner core of the Earth may spin faster than the 
rest of the planet so that the above exact solution does not hold for 
the fluid motion of the terrestrial core. This conclusion implies that 
the temperature distribution can not be a radial distribution due to 
the fact that Eq. (6) is, in general, coupled with the velocity field. To 
enforce this statement we shall show that in general a rotating fluid 
sphere, as previously stated, can not be described by the <f> component 
of velocity. Let us then assume that a sphere starts to rotate with 
Vfp{t) such that ^^(0) = and vg = Vr = 0. We also assume that for 
symmetry all physical variables do not depend on (j). Writing only the 
hydrodynamic part of the set of equations (1) - (4), i.e. setting B = 0, 
we obtain 




(11) 
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(12) 



^ idP 

r p dr 

^ v^. (13) 



dt p 
Combining Eqs. (11) and (12) we obtain 

The above equation is satisfied by a function tJ<^(r, 9) of the form v^{r, 9) = 
vj){rsin9). Using this result we may rewrite Eq. (13) as 



dt p \ dx"^ X dx X 



Vci>{t,x) = -\^ + -—~-^\v^{t,x), x = rsui9. (15) 



The general solution of the above equation, via Laplace transform and 
with the condition Vfi){Q,x) = 0, is 



v^{s,x) = ah\-^\+hKi\^\ (16) 




where, by definition, ?)^(s,a;) = exp[— st]t;0(t, x)(it and D = rj/p. 
Assuming a finite solution for r — t- and consequently a; ^ 0, then h = 
0. It is evident by the form of Eq. (16) that the boundary condition 
for the velocity, i.e. ^^(s, i?sin^) = ali ^i?sin where R is the 

sphere radius, can be satisfied only by a restricted class of functions. 
One could consider adding another component, for example the ra- 
dial component Vr-, to describe matter falling in the center. But the 
independency of the dynamic by the angular variable (f) implies that 
the 9 component of the velocity vg also has to be considered. This 
fact comes directly from the continuity equation. Let us assume that 
there is another component of velocity, the radial component Vr- The 
continuity equation is 

^2^y-r)=^. (17) 

The solution Vr = f{9)/r^ diverges at the origin and can not vanish 

on the surface of the sphere r = R. To avoid this inconsistency we 
are forced to add vg to the flow. From this we can infer that from the 
early stages of Earth's formation to the present, the velocity of the 
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fluid could not be described by only one component of the velocity 
vector. This analytical conclusion is in agreement with the numerical 
works presented in the references. 

4 Thermal generation of magnetic field 

Looking at Eq. (3), or Eq. (8), we notice that is a diffusive-like equa- 
tion. For the time evolution of B it is important to give an estimation 
of its initial value. It is accepted that the Earth's core is made mainly 
of iron with a solid inner core the size 10^ km and an outer core of 
liquid about 2 x 10^ km thick [7] . The temperature distribution of the 
core is not uniform and it ranges from approximatively 10^ °K at the 
very center to 10^ °K at the surface of the outer core. The non uni- 
form temperature can generate a contribution to the electrical current 
that is proportional to the gradient of the temperature, known as the 
Peltier— Seebeck effect, so that we can write the total current as [13] 

J = (7[E + vxB - a(T) VT] . (18) 

Note that Eq. (3) docs not change if a;(T) VT can be written as a gra- 
dient of a function. To evaluate the coefficient a, we have to consider 
the fact that the density of either the solid inner core, or the fluid 
outer core, is such that the electrons can be considered a degenerate 
Fermi's gas. Indeed, according to the authors of Ref. [14] the core 
density is of the order of 10^ Kgm^^. This implies that the Fermi 
energy of the electrons of the Earth's core is 

Ii2 / AT\2/3 

eF = — (Stt^—) f« 2 X 10-^^ Joule (19) 

corresponding to a Fermi temperature Tp ~ 1.4 x 10^ °K. This is at 
least one order of magnitude higher than the Earth's core temper- 
ature therefore justifying the degenerate Fermi's gas approximation. 
Quantum calculations show that [13] 

I a(r) h (20) 

where kB is the Boltzmann constant and e is the electron charge. Let 
us consider a very simplified model of the early stage's of Earth's for- 
mation. Models suggest that Earth's core was completely molten [15]. 
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Since the temperature change of the core has a time scale on the or- 
der of Earth's age [16] we infer that whatever is the contribution from 
the thermal term in Eq. (18) this contribution still holds today with 
the same order of magnitude. A widely accepted estimation of the 
core temperature is approximatively 8 x 10^ °K for the inner core 
and 4x 10^ °K for the outer core (see for example [14, 17, 18]). Let us 
consider the contribution to the magnetic field due to thermal term. 
Prom Maxwell's equation we obtain 

V X B = iiaJ = -^oo-a(T)VT. 

In general, the temperature distribution in time and space is coupled 
with the velocity field, Eq. (6), and as shown in the previous section, 
all components of the velocity are present so that the temperature 
distribution can not be only radial. Note also that we use fiQ as the 
value of the magnetic permeability since at such temperature we as- 
sume that there is no magnetization. We can deduce the field scale 
via the relation 

-«/xoaa(r)— (21) 

where AT is the difference in temperature, and R is the length scale 
of the system. We obtain the scale strength of the thermal magnetic 
field 

Bt = iiQG^^^kB^T. (22) 

eep 

Using the value of temperature Tc « 8 x 10^ °K, AT ps 4 x 10^ °K, 
and cr 10^ S m~^ [19] we obtain the numerical value for the strength 
of the core of Earth 's magnetic field 



Bt ~ 0.0024 tesla = 24 gauss. (23) 

The estimated strength of the core of Earth 's magnetic field is ap- 
proximatively Best ^ 25 gauss [20] which is very close to the analytical 
value given by Eq. (22). We note that selecting different values for the 
temperature, according to the different models present in the litera- 
ture, the value of the field would of course change consequently, but 
the scale magnitude remains of the order of tens of gauss. The ther- 
mal current gives a strong (if not total) contribution to the magnetic 
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field. It is worthy to stress that in principle the only phenomcnological 
parameter, i.e. the conductivity cr, could be evaluated using quantum 
mechanics [21] so that we can conclude that the field given by the 
expression (22) may be written in terms of fundamental constant and 
physical parameters of the system such as density N/V and temper- 
ature T. Note that dependence on the radius R of the fluid region 
is implicit in the dependence of the temperature on R. This is the 
main reason why we kept explicit the temperature difference AT in 
Eq. (22). Once we clarify this, we can rewrite more concisely Eq. (22) 
as 



where T can be taken, for example, as the average temperature of the 
fluid region. We conclude that Eq. (24) represents the scale of the 
strength of the magnetic field of celestial bodies with an Earth-like 
physical condition for the core, from a statistical point of view. 



5 Jupiter's magnetic field estimation 



In principle we can apply the ideas of the previous section to other 
celestial bodies, particularly in our solar system. The main difficulty 
with this is the scarcity of information about the physical internal 
condition of other planets, although we can make some general con- 
siderations. For example. Mercury and Mars are quite smaller than 
Earth. This fact surely contributes to a faster cooling of their interiors 
so we can expect that the cores of these planets are no longer in the 
fluid state. In fact the two planets have a very weak magnetic field. 

Venus does not have a magnetic field [24] and there are several 
possible explanations for this. Venus is a planet very similar to Earth 
in dimension but it docs not exhibit volcanic activity, and this could 
imply a cold core. Venus has a very slow rotational motion compared 
to Earth, and a rotational motion is considered to play a crucial role for 
the terrestrial magnetic field. Also we should consider the possibility 
that Venus could be in a reversal phase. 

Jupiter is a good candidate to test our model. Even though little 
is known about the planet, its internal structure has been modeled by 
several authors (see for example [25, 27, 26]) and the physical informa- 




(24) 
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tion is enough to allow a rough estimation of the scale of its magnetic 
field using Eq. (22). The electrical conductivity is (t f« 10^ Sm"^ [26], 
its temperature ranges from T 2 x 10^°K for the core boundary to 
T ~ 10^°K for the metallic hydrogen boundary, and the estimated 
density of the metallic hydrog 4 X 10^ Kgm"^ [25, 27]. Con- 

sequently the Fermi's energy takes the value ep ^ lO"-*^^ Joule corre- 
sponding to a Fermi temperature Tp 7x 10^ °K so that we can apply 
the Fermi statistic for the electrons in the metallic region. Plugging 
these values into Eq. (22) we obtain for the magnetic field of Jupiter 
an estimation of its strength in the metallic hydrogen region Bj k, 30 
gauss. Taking into account that this region extends for a fraction, 
ranging from 0.7 to 0.78 times the Jupiter radius, we obtain for the 
surface value of the magnetic field a range from Bjs « Bj(0.7f « 10 
gauss to Bjs ~ 5j(0.78)^ « 14 gauss. This is in agreement with 
the scale of the observed values [28]. As said for Earth, according to 
several models, we can change the values of the parameters, but the 
strength of the field would still be of the order of the observed field. 

6 Conclusion 

We provided an analytical estimation of the magnetic field scale of 
planets with physical core conditions similar to that of Earth from 
a statistical point of view. The magnetic field strength was evalu- 
ated directly from the physical parameters of the center of the planet, 
considering density, temperature and core's size. We showed that 
an electrical current generated by a thermal gradient can exist in a 
rotating fluid sphere and can give an important contribution to the 
magnetic field. Our conjecture was supported by estimating the mag- 
netic field strengths of Earth and Jupiter that were in agreement with 
the observed magnetic field intensity of the two planets. 
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